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ABSTRACT 


The two-phonon processes in slow neutron scattering against single crystals are discussed. 
It is shown that in the harmonic approximation the energy distribution of neutrons scattered 
coherently in a certain direction in some cases has logarithmic singularities for the energies ob- 
tained in the one-phonon processes. Further there appear singularities of the same kind as Van 
Hove found earlier for the incoherent one-phonon scattering. To illustrate the logarithmic 
singularities more explicit results are given for an ideal Debye crystal. A formula for the incoherent 
two-phonon scattering is also given, and this is applied on KCl. 


I. Introduction 


In recent years several successful experiments have been performed to determine 
the eigen vibrations of crystals, using scattering of slow neutrons [1-5]. One there 
uses the fact that in processes where one vibrational quantum participates the 
neutrons are scattered with a certain number of discrete energies, which are well 
separated and can thus be measured. Such experiments will certainly be a valuable 
complement to those on diffuse X-ray scattering [6-9] in providing more complete 
information about the crystal vibrations. So far much stronger intensities are obtained 
using X-rays, but owing to the fact that the vibrational frequencies are determined 
from the measured scattering intensities, corrected for Compton scattering, higher 
phonon scattering etc., the interpretation of the experimental results is not as direct 
as for neutrons. When using neutrons only the scattering energies are measured, and 
for a well monochromized incident beam the intensities mostly have sharp peaks for 
energies, corresponding to one-phonon scattering. Owing to incomplete monochroma- 
tization and collimation the peaks always have finite widths.1 We cannot @ priori 
assume that especially the coherent two-phonon scattering does not display more or 
less pronounced intensity maxima and minima.’ If such exist they may in some cases 
confuse the interpretation of the experimental results. It should also be noted that 
the incoherent one-phonon scattering may have a similar confusing effect.® 


1 Van Hove has shown that even under ideal experimental conditions the one-phonon scattering 
surfaces have finite widths, owing to the anharmonicity of the crystal vibrations. 

2 Brockhouse and Stewart [2] claimed to have found indications of such maxima in their 
experiments. Carter, Palevsky and Hughes [4] also got some peaks that could not be interpreted 
as one-phonon peaks, but as Palevsky informed the author these results were very uncertain. 

3 See G. Placzek and L. Van Hove, Phys. Rev. 93, 1207 (1954). 
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So far very few discussions of the corrections to the one-phonon scattering due 
to higher phonon processes have been published. General expressions for their in- 
tensities are given by Fréman [10], Glauber [11] and Van Hove [12].? It is, however, 
hard to draw any further conclusions from their formulas without a more detailed 
analysis of these expressions. Several authors [13-16] have discussed the effect of 
the higher phonon processes on the total cross section, but here we are interested in 
the differential scattering cross section, indicating the energy distribution of the 
scattered neutrons. Similar problems arise in connection with X-ray diffraction, but 
owing to the fact that different cross sections are considered in the two cases the results 
on X-rays cannot be applied to neutrons [6-9]. 

In the one-phonon processes the energies of the scattered neutrons are correlated 
to the eigen frequencies of the crystal by the following energy-momentum rela- 
tion [17]: 

2 2m 
=k —2 =" wy (9), Ks, 
k=k)+q-t. 


ky, k and q are the wave vectors of the incident and scattered neutron and the wave 
vector of the eigen vibration, participating in the scattering. t is an arbitrary re- 
ciprocal lattice vector, defined in such a way that the product of this and a direct 
lattice vector is a multiple of 27. w;(q) denotes a certain branch of the eigen frequencies 
and m is the neutron mass. ¢ can take the values +1, corresponding to phonon absorp- 
tion and phonon emission respectively. 

In the following we will discuss the two-phonon processes to determine the charac- 
teristic features of the scattering. Our main aim is to ascertain whether the cross 
section shows any singularities and to get an idea how sensitive it is to the details of 
the crystal vibrations. 


II. The high phonon processes 


For the low order phonon processes we know that the coherent and incoherent 
scattering cross sections vary in very different ways with the scattering angle and 
energy, the coherent ones showing stronger singularities. Concerning high phonon 
processes it can, however, be shown that the two kinds of cross section are approxi- 
mately proportional. The intensity from the n-phonon processes is approximately 
given by? 

Fo" yak e-ev OW) 
dQdw ko n! Gy (22n)* 


bn(w—a), n>. (2) 


N is the number of nuclei in the crystal. a should be set equal to the coherent or 
incoherent scattering length respectively, and e-*” denotes the ordinary Debye- 
Waller factor. iw, and fiw are the energies of the incident and scattered neutron 


? Some reults are given by P. A. Egelstaff in Harwell Reports AERE-1164. Some remarks 
may also be found in: G. Plazcek and L. Van Hove, loc. cit. 


* The formal treatment goes back to ideas found in A. Akhieser and I. Pomeranchuk, J. Exptl. 


Theoret. Phys. U.S.S.R., 17, 770 (1947), though these authors give explicit formulas only in the 
incoherent approximation. 


% The proof will be published in a later paper. 
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and bn (w— @ ) is a normalized Gaussian function with the width a, Vn n 


dy (2 an)? 
and the maximum value at w—w)= —nm, o, and m being temperature depend- 
ent constants independent of n. From (2) we thus conclude that when 2W <1 the 
intensities from the higher phonon terms decrease rapidly with increasing order. We 
can also assume that the higher order scattering varies much more smoothly than for 
instance the second order scattering. In typical one-phonon measurements using 
cold neutrons even the two-phonon intensities are rather weak, and we may there- 
fore neglect all higher phonon terms. ‘ 


IIIf. The coherent two-phonon processes 


For brevity a Bravais’ lattice only will be considered. As a starting point we may 
take the following expression for the coherent two-phonon cross section per unit cell, 
given by Froéman [10]:! 


2 : h (x- C;.)? 
Soda = 7 ko oe rairs 2NM wy 


th 3) + th v) te’ ; 
Ee: OE PONTE 6, (e+4+4') 8 (ay tea, +6’ de). (3) 
Some of the notations have been explained above and the new ones denote the 
following: 


€, =polarization vector of an eigen vibration, characterized by the double index 
s=(q,7).? 

M = mass of a nucleus. 

V =volume of the elementary cell of the crystal. 

B =h/2k,T, where kz = Boltzmann’s constant, 7’ = temperature. 

x =k, —k. 


6(...) is an ordinary Dirac function and dp (x =2 6 (x — 7) denotes a periodic Dirac 


function singular at the reciprocal lattice points. - double summation should be per- 
formed over the eigen states s = (q,7) and s’ = (q’, 7’) and a summation should also 
be made over the two values -++1 for ¢ and ¢’, 

The Dirac functions express the conservation of energy and momentum, analogous 
to the case of the one-phonon scattering (see Fig. 1). 


gar ,2m ? 
os w;(q) -€ ct, wy’ (q if 


k=kj+ q+ q'—*. 


aps (4) 


We note that two phonons are now involved. 


1 The notations used here are not exactly the same as in Fréman’s paper. 
* (, is normalized according to |¢, | = 


bo 
— 
~I 
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ko 


Fig. 1. 


Before any conclusions are drawn from the above intensity formula, the latter will 
be brought into a more convenient form. To save space we introduce the notation 


(ey ‘ C;) 


2 
fi (4, €) = - [coth (6 a.) + e}; (5) 


. Ws 


ex then means the unit vector along x. 
The summations over q and q’ are transformed into integrations over the wave 


vectors as 
NV 
saath 
=> Gap } i 


(cell) 


and new integration variables q, and q, are introduced as follows: 


qi = 4, G=qtq’. 


The 6p-function only contributes when (x + q,) reaches a reciprocal lattice point and 
after the q,-integration we can write 


AU eee Eke RPL a Kt oY 
lear) (2 2)? 2 ~ 2 * 


rt” 8,6" 


PRO Viet! Meee 


x WACT é) fy(tT—*— qy, €’) 0 [@ — Wg + € @;(q,) + €' oy (tT — 2% —q,)] dq. 
Ve 


> denotes the summation over the eight possible lattice points above, and the 
Tv 


volumes of integration are determined from the fact that both qi and (q, — q,) 
should be contained in a reciprocal unit cell. The volumes V~ do not overlap and taken 
together they precisely fill a unit cell. Using the periodictiy of f,(q,¢) and also noting 
that f;(—4q,¢) =f,;(q,¢) we can simplify our expressions as follows:1 


> f-dq= [tan €) frde+ dy, €) d[@— oy + £05 (q) +6" oy (e+ 4,)] dq. 


Tt Ve (cell) 


1 This latter transformation is used only to get a more attractive form. 
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The 6-function stipulates that only those q,-vectors should be considered which lie 
on certain surfaces in the unit cell, defined by: 


Let us call such a surface a A-surface. We note that we get one A-surface for each 
combination of the branches 77’ and also for each combination of ee’. Further the 
A-surfaces vary with x, a fact which makes the numerical calculations very tedious 
when the intensities for many x-vectors are required. 

The above 6-function can be eliminated by transforming to a surface integral. 
Thus denoting a A-surface by S, and with gradA meaning the gradient with respect to 
q we may write” ” 


a 2), ee k e-2wl (55 


2 
ais coh 
d @) dw —_ Acon ko 9 2 a) G3 ( if 


ws 2 72)° 45 | grad A| 


GS" ( )= us >: | fi (q, é) fy (x+4q, e') ag. 
ne, 


When determining the possible singularities of the cross section we must pay 
special attention to the following critical points of the integrals: 


(1) q=0, q=t—x, in which points f;(q,¢) or f;(% + q,¢’) are infinite. 
(2) zero-points of grad A. 


The intensity formula shows some characteristic similarities to the corresponding 
expression for the incoherent one-phonon scattering [18].° 


2 (1) 2 
d Oincoh 2 k _opwhx 


incoh pi 
dQ a ra maul ) um” (w Wo), 
8 
Gincon (co) = {) feoth (Bw) —1], |w|<1, (8) 
2a | 
= 0. Jol>1. J 


f(w) means the vibrational frequency function with the neutron “‘frequency” 
inserted instead of the vibrational frequency, and f(w) is defined for negative fre- 
quencies as follows: 


{(—@) =f(o) 
1 

and normalized to unity* / f(w)dw=1. 
0 


1 In S; and grad A we have for brevity omitted the indices jj’ ee’. 

2 It must be noted that for cold neutrons the A-surfaces corresponding to emission of one or 
two phonons do not always exist, which implies that the corresponding integrals vanish. 

3 The formula is exactly valid only for simple crystals with cubic symmetry. See also G. 
Placzek and L. Van Hove, loc. cit. 

4 The frequencies are measured in units of the maximum vibrational frequency. 
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To make the resemblance more obvious Gj"°°"(w — ) is written in the following 
altered form: 


i us 9; (4; €) 
incoh sar = 9 
CY { (Oreo) 3 (2 2)° 2 | | grad A | oe: (9) 
Sa 
where g; (q, €) = e [coth (6 w,;) + é]. (9’) 


The integration should now be performed on the iso-frequency surfaces 
Aj,e = — Wy + €w;(q) = 0. (9’’) 


We know that the incoherent one-phonon cross section varies continuously and 
that singularities of the following type appear [19]: 


a? ofn 
incoh _ 4 

qOala =const|w—a,|', o>, 

0, Wo < We. (9"””) 


These are due to the vanishing of grad 4. Moreover for some crystals the intensity 
shows sharp maxima and minima.! 

Owing to the resemblance of (7) and (9) we might draw similar conclusions for the 
coherent two-phonon cross section. To ascertain the existence or non-existence of 
singularities in the two-phonon intensity let us consider the critical points in more 
detail. 


IV. The critical points q=0 and q=1—x 


From (6) we see that when a A-surface passes through either q =0 or q=t—x, 
where t is chosen so that q lies within the reciprocal unit cell, the momentum vector 
and the energy of the scattered neutron are correlated in accordance with the one- 
phonon scattering conditions. Let us consider a certain branch of the one-phonon 
scattering surfaces, which thus is determined from 


W@ — Wot €0;, (x) = 0 
for definite ¢, and j,. When this condition is fulfilled the A-surfaces indicated by 
(77, € €') = (Jy, € €1) 
for arbitrary j and e go through q=0. Correspondingly the A-surfaces with the 
indices 
(17, €€') = (919, ee’) 


for arbitrary j’ and e«’ go through the point q=T—x. 
The two critical points are completely equivalent and therefore we need only 
consider one of them, say q = 0. To determine how a corresponding integral included 


* See for instance the frequency spectra drawn by Jakobsen [8] for Cu and by Walker [9] and 
Squires [20] for Al. 
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in (7) varies when x approaches a point x, on the one-phonon scattering surface, x 
is supposed near x, outside the surface (Fig. 2). In that case the A-surface is approxi- 
mately determined from the equation below, if we keep near q = 0. 


A=ec,(q)|q|+ 6’ grad w;(x,)-q —u =0, (10) 
where c¢,(q) denotes the sound velocity. uw is defined by 
| [Vv are. grad WM; (*.)]- (* =x,); (10’) 


where v denotes the velocity vector of the scattered neutron. 


[ioe 2a 


Suppose that grad 4+ 0 everywhere on the A-surface. Otherwise we would have 
an even stronger singularity. Except for a strictly isotropic crystal c; depends on 
the direction of q, and we may have the following two cases: 


1. |grad w, (%,) |< inf c;(q): 
The integration surface is contained within a sphere with radius r ~u, and we have 
mainly to evaluate the integral 
dS 


Lo; (q)] 


A 


which tends to a finite value, when j.—>0. We can also show that the intensity tends 
to the same value when approaching the one-phonon surface from both sides. This 
case will surely occur, for grad w,(x%,) actually vanishes at certain points, and for 7 
denoting the longitudinal and j’ a transversal branch c, will ordinarily be larger than 
|grad w;| that is when the one-phonon surface corresponds to transversal vibra- 
tions. 


2. |grad w; (,)| > sup ¢,(q): 
The J-surface has the shape indicated in Fig. 3 and we may proceed as follows, 
excluding all constant factors and only considering the order of magnitude:t 


(00) 
dS dS | rdr =| 
ge -(*%+q, e) ——— ~ —_ w | —|. 
ic ) fr ( 1 'Terad al la Pe u 
S) Sy 0 


1 y is inserted in the figure as the radius, using cylindrical coordinates. 
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Fig. 3. 


The cross section will thus be infinite when x approaches the scattering surface 
according to the formula 


& 62 
ee —In |[v—e’ grad a; (*,)] + (% — %-) | - (11) 


we may also have the same type of singularity, but the exact condition for this is 
rather complicated and will not be given here. This case may especially happen 
when j’ refers to the longitudinal and j to a transversal branch or in other words if 
the one-phonon surface refers to longitudinal vibrations. Experimentally the 
singularity would appear as a broadening of the one-phonon scattering surface. 


Considering the above derivation in more detail we may conclude that for scattering 
points approaching a reciprocal lattice point t along a one-phonon surface the inten- 
sity increases as 1/|% — |? besides the logarithmic increase above. 


V. The critical points grad 2 = 0 


The zeros of grad A are critical points of our integrals only when they lie on the 
integration surfaces, and the conditions for this (A = grad A = 0) may be written in 
the following form: ; 

EM; (q) + &’wy(x + q) +m —@, =0, 


12 
€ grada@;(q) +’ grad w,(x + q) =0. a 


These equations determine a certain surface Sy, with the following property: 


For a x=x, on Sy there exists at least one critical point q, of the corresponding 
integral. 


A x% near Sy is now considered and the contribution of the integral from the region 
near a critical point is estimated. Supposing the critical point to be analytic the 
A-surface may be approximated by a surface of second degree, when keeping near q,. 
The equation of the surface can be written as 


dy Ef + de 3 + a3 63 =i, (13) 
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where § = q — q, and the coordinate axis have been chosen along the principal axis. 
@,, Az, 43 are positive or negative constants and y, is defined by 


by — Lv at grad (Ke as q-)] *(* AE (14) 


It is now a straightforward matter to evaluate the integral (see [21]), and this is 
shown to have a singularity of the following type, when crossing the Sy-surface:1 


ec const |[v—e’ grad w; (x%,+q.)]*(*%—%,)|* on one side 
=0 on the other side of Sy,. (15) 


[vy — «’ grad w,(%, + q,)] is normally directed to Sy. The above singularity corresponds 
exactly to that for the vibrational frequency function. As the same kind of discussion 
may be found in other papers [19, 21] on the frequency function, we refer to these 
for a more complete treatment, including also non-analytic critical points. 


VI. The singularities for an ideal Debye crystal 


Though the scattering cross section seems to be rather sensitive to the detailed 
assumptions on the crystal vibrations, it may be illuminating to apply our results 
to an ideal Debye crystal. We then have 


1 =W, =W3 =@(q) =c|q| (c¢ =sound velocity), 


if q is restricted to the reciprocal unit cell,? and the A-surfaces defined in (6) are the 
same for all values of (j,7’). We can thus sum over these indices in formula (7) and 
write the result in the following form: 


I 1 
oh i (q, €) fr (q's &) = eP (@o- ®) 5 
79" By 


c (q) sinh [8 w (q)] 2 @ (q’) sinh [f @ (q’)] 


with q’=x + q. 
For sufficiently high temperatures sinh[fm(q)] can be approximated by 1/fa(q) 
which further simplifies our expression. We get 


a V oat 1 ds 
he} (2x) ce J laf? |xta—tP [grad A| 
S; 


GS (x4) = EP (Mo — ®) ( (16) 


where t is chosen so that (x + q—7) always lies within the unit cell. The /-surfaces 


obtained from 
Wy — @ 


e|q|+e’|x+q-—7|= (17) 


are either ellipsoids or hyperboloids or parts of these with their foci at q =0 and 
= —x depending on the position of the scattering point. We will especially deter- 


1 The continuous part of the cross section has been excluded in (15). 
2 In the Debye approximation the unit cell is spherical. 
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mine the scattering near the one-phonon scattering surfaces. Owing to the symmetry 
we need only consider one of these, for instance that corresponding to phonon emis- 
sion (@<«@,). When we keep on one side of the surface w)—w —c |e—x|>0, 
and according to (17) only the integral, indicated by (¢,¢’) = (1,1), contributes to 
the intensity. On the other side of the surface wy — @ — c|t—x|<0, and in that 
case the non-vanishing integrals correspond to (¢,e’) =(1, — 1) and (—1,1). Fora x 
very near the scattering surface only those parts of the A-surfaces passing near q = 0 
or ¢=t—x! need be considered, and the equations for them may be written in 
the form 


xv yte 
ee aI 18 
at B (18) 
2 
: 2_ (Vo @ 
with i= ays ); 
: ! e (18’) 
2 Wo @ Ti 
: a + rYRM 


The x-axis has then been directed along the line between the two foci and the origin 
(x =0) is placed midway between them. The following treatment will be separated 
into two cases as b7>0. 


‘bj >0, corresponding to scattering points on one side of the one-phonon surface 
(Fig. 4a): 
We may proceed as follows when evaluating the integral corresponding to (¢,é’) = 
(1,1). The others vanish. 


a, 


ff 1 Te [ % ae 
J laP Tera sP fered 2)? J [uP Tet a—eP [(erad 2), 
A = % ' 


a, 


ie 22 dx 
oo J [allerq—s| 


__2ne fee l Ja 
(o—o) J [fa] |xt+q—|]°" 


-—a, 


ane ‘da 


(@) — oy la] 


r= (y? +2")? and (grad J), denotes the component of grad 4 perpendicular to the 
x-axis. We note that for g running along the A-surface 


|x—t| 2 
= iP —_=_ — 
|q| a( + Tres 


1 t should then be chosen to restrict (T—%) to the unit cell. 
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——_ 


—__—— 


Fig. 4a, b. The integration surfaces. 


and the final integration thus gives the following result: 


| --as- ae In 
|x—x| (@y — w) 
Sy 


Wy- ot+e|t—x|| (19) 


Wy— o—c|t—x| 


bj <0, corresponding to scattering points on the other side of the one-phonon surface 
(Fig. 4b): 
We ae in a similar way as above. Only the integrals corresponding to (¢,¢’) = 
(—1,1), (1, —1) do not vanish. We now get 


(00) 
2 dS 4c (el i “| a 
x 
airires <P [grad A] (o—o)°. |xt+q—t| |q| 
and further 


|c—x| x [ce —x| x 
lal=a( fee tl)) eta [=o (5 2-1). 


1 


Extending the integration to infinity we obtain the same result as above. 

Our formula will also be valid near the scattering surface for phonon emission if 
we only substitute |w)—@| for (wy —@). 

The differential scattering cross section can now be written in the following form, 
valid in the Debye approximation and near a one-phonon scattering surface: ** 


d? o®), * k =i h2 x 67 2 
coh _ a, oil errr erat bad : 
oe = Foon exp[—2W+ f(a —o)l 5 \2 M kz 0 =) a (20) 
] w?, | wy — w| + @ (x) 
iron J , 
(x) 12m w(x) |@—o| Perens 


Here 6 = Debye temperature; w,, = k,9/h =maximum vibrational frequency. 


1 Note that for cold neutrons the one-phonon surfaces corresponding to phonon emission 
disappear. 

2 Egelstaff (loc. cit.) considered this case only and his result essentially corresponds to formula 
(20), though he did not give this explicit expression. 


225 


A. SJOLANDER, T'wo-phonon processes in neutron diffraction 


VII. The incoherent two-phonon scattering 


To complete the discussion on the two-phonon scattering we also give the cross 
section for the incoherent scattering. Most crystals certainly scatter neutrons coher- 
ently much more strongly than incoherently, but for some crystals, as for ete 
vanadium and crystals containing hydrogen, the incoherent scattering dominates. 


2.0 1,0 0 1,0 2.0 
Ge ( w) 


Fig. 5. To get the normalization propo*cd the values obtained from the figure must be multi- 
plied by the factor 0.954. 


The scattering cross section may be expressed in the vibrational frequency function 
f(m), and for convenience we introduce the following notation:? 


f(o) 
2m (0) 
==(), |w|>1, (21) 


1 
where y(0) is defined from? J g(w)do=1. 


g(w)= [coth (Bw)—1], |w|<1, 


=! 


1 Recent experiments on vanadium have been made by Brockhouse [22], Carter, Hughes and 
Palevsky [23], and Stewart and Brockhouse [24]. Results from similar experiments on ZrH are 
found in refs. [26, 27]. 

2 See formula (8) and the text below. The frequencies are measured in units of the maximum 
vibrational frequency. Compare also A. Sjélander, Ark. Fys. 13, 199 (1958). 

° Note that 2W = (hx2/2M) y(0) 
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The cross section can then be written in the following form:! 


477A is the total incoherent cross section for a single nucleus. In the integration over 
«, most of the details in the frequency function disappear and thus the incoherent 
two-phonon scattering ought to vary more smoothly than the coherent one. G°™ 
has been computed for KCl at room temperature, using numerical data on the eigen 
vibrations from Tenerz [25], and the result is given in Fig. 5.? 
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